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Abstract 



A theorem about asymptotic estimation of multiple integral of a special type 
is proved for the case when the integrand peaks at the integration domain bound, 
but not at a point of extremum. Using this theorem the asymptotic expansion of 
the electromagnetic deuteron form factors at high momentum transfers is obtained 
in the framework of two-nucleon model in both nonrelativistic and relativistic im- 
pulse approximations. It is found that relativistic effects slow down the decrease of 
deuteron form factors and result in agreement between the relativistic asymptotics 
and experimental data at high momentum transfers. 



PACS numbers: ll.lO.Jj, 13.40Gp, 13.75.Cs 
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1 Introduction 



Recent advances in experimental investigations of hadron structure arouse the interest 
in the theoretical study of the hadron electromagnetic form factors at high momentum 
transfers (see e.g. review [i| and references therein). In this connection the JLab program 
of investigations on elastic electron-deuteron scattering experiments at ^ 10 (GeV/ c)^ 
(Q^ = —q^ , q is transferred momentum) [2] attracts considerable attention. There exists 
a hope that these JLab experiments will help to determine the limits of application for 
the two-nucleon model and to clarify the interplay between nucleon-nucleon and quark 
approaches to the deuteron. 

Using the asymptotic expansion presented in this paper we show [3] that the momen- 
tum transfer region in the JLab experiments is asymptotical for the deuteron considered 
as a nucleon-nucleon system. That is why the study of the electromagnetic deuteron form 
factors is interesting at — >• oo. 

The present work is devoted to the theoretical investigation of deuteron form factors 
asymptotic behavior at high momentum transfer in the framework of the two-nucleon 
model. The form factors asymptotics is studied in accordance with the next points. 

1. As a rule, calculation of the form factors asymptotic behavior in the relativistic 
approaches reduces to the asymptotic estimation of n-tuple integrals. In the relativistic 
approach used in our work the deuteron form factors are expressed in terms of double 
integrals where integrands peak at the integration domain bound, and corresponding point 
is not a point of extremum. In this connection the theorem defining asymptotic expansion 
of n-tuple integrals with such integrand is proven in our paper. 

2. In general, high momentum transfers require relativistic consideration. However we 
begin the consideration of the asymptotic estimation of electromagnetic deuteron form fac- 
tors with the nonrelativistic case and nonrelativistic impulse approximation at oo 
. This is due to the facts that, at first, the nonrelativistic calculation is a less complicated 
and, second, this calculation is important for establishment of the role of the relativistic 
effects. 

3. The asymptotic expansion of the relativistic deuteron form factors is calculated in 
the relativistic invariant impulse approximation in a variant of instant form of Poincare- 
invariant quantum mechanics (PIQM) developed in our papers previously [H O [6l [TJ [8] . 
The relativistic calculations are performed by analogy with nonrelativistic case. It is 
shown that relativistic effects essentially slow down the asymptotical decrease of the form 
factors. 

4. It is found that obtained in the framework of the two-nucleon model relativistic 
asymptotics coincides with the experimental data. 

This paper is organized as follows. Section 2 is devoted to the proving of the central 
for this work theorem defining asymptotics of multiple integrals of some special type. In 
Sec. 3 a brief review of the formulas for the deuteron form factors in the nonrelativistic 
and relativistic invariant impulse approximation is given. The deuteron form factors 
asymptotics is calculated in nonrelativistic and relativistic impulse approximation with 
the help of the proven theorem in Sec. 4. In Sec. 5 asymptotics of the form factors is 
obtained for the deuteron wave functions in the conventional representation as a discrete 
superposition of Yukawa-type terms [9J. Sec. 6 contains the conclusions of this paper. 
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2 Theorem on the asymptotic expansion of some mul- 
tiple integrals in the case when the maximal value 
of the integrand belongs to region boundary 

In the following we will consider integrals of the kind: 

F{X) = f f{X,x)exp[S{X,x)]dx , (1) 



where f2 is a domain in R", ), A is a large positive parameter. We will use 

following definitions: dQ is a bound of the domain Q, [Q] = QU dQ, the bound dQ G C°° 
if in the vicinity of any point x° G dQ it can be specified by equation Xj = (p{x'), x' G 
?7', x' = (xi, ...,Xj-i,Xj+i, ...,Xn), U' is a neighborhood of a point x °, and the function 
^{x') G C~ in U'. 

Let us prove now a lemma and a theorem on the asymptotic estimation of integrals in 
Eq.©. 

Lemma. Let S{X,x) be a smooth function in the [Q], f{X,x) be a continuous 
function in the [Q], and M(A) G C^, 

M(A) = sup S{X,x) < oo , 

at some Aq > the integral ([T]) be absolutely convergent: 

|/(Ao, x) \ exp[S'(Ao, x)]dx < oo , 



n 



and the following conditions be fulfilled at X > Xq: 

dS{X,x) dMjX) 

dX - dX ' ^ ' 

\f{X,x)\<C^\f{Xo,x)\. (3) 
Then at X > Xq the following estimation is valid: 

\F{X)\<C2e^'^^K (4) 

Proof. At A > Ao the following estimations are true: 

|F(A)| < e^^W /■e5(Ao,x)-M(Ao)gS(A,x)-5(Ao,.)-M(A)+M(Ao)|^(;^^^^|^^ < 



< gAf(A)-Af{Ao) J ^S{Xo,x)+S{X,x)-S(Xo,x)~M{X)+M(Xo)^j^(^^^^^^^^^_ 

n 

From conditions Q,© we obtain the inequality: 

|i^(A)| < Cie'n^)-^^(^o) f e^(-^<^'^^\f{Xo,x)\dx < Cse^W . 



n 



Thus the statement (jlj) of the lemma is proved. 

Later we will consider function ^(A, x) described in lemma which has the maximal 
value in the point x° G dQ, and ^(A, x), dQ G in the vicinity of x^. This point is not 
the point of extremum, that means the validity of the following conditions: 



dS{X,x'^) 
dn 



7^0 



and matrix of coefficients B: 



d^S{X,x^ 



n-l 



B 



(5) 



(6) 



gives the negative determined quadratic form. In Eqs. d/dn specifies the internal 

normal derivative n to the dVL, and ^i, is an orthonormal basis in the tangential 

to the dVL plane T dVL^o at the x'^ point. 

For convenience let us choose in the vicinity of point x^ a frame y = [yi, ...,yn), so 
that x° is the origin of coordinate and the internal normal to dQ coincides with the last 
basis vector of the new coordinate system. Functions /, 5* in these coordinates we denote 
as /*, S*, and U* is an image of U (that is an image of a half- vicinity of the point x^). 
The equation for dU* in the vicinity of the point y = can be written in the following 
way: 

Vn = (p{y'), y e U' , y = {yi, z/„-i), 

with U' is a vicinity of the point y' = 0, (p{y') G C°°([/'), and at y' 0, ip{y') = 0{\ y' p). 
Theorem. Let the following conditions he fulfilled: 
1°. /,5GC([fi]). 

2°. ^(A,^^) is maximal value of function S'(A,x) , x^ G dVL, and x^ is not point of 
extremum. 

3°. /, S*, dVt G C°° in the vicinity of the point x^ . 

4°. The Taylor expansion of functions S* and f* in the vicinity of point x^ satisfy the 
following relations: 

r(A,y) = r(A,o)[i + o(i)], (7) 

S*iX,y',^iy'))-S*iX,0) = ^{Ay', y') + 0{\ y' \') , (8) 



the matrix A = 
Then at X 



n-l 



angle brackets denote bilinear form: {x, y) = Xiyi + X2y2' 
oo the following asymptotic expansion is valid: 



F(A)~exp[5(A,x°)]^a,(A) 

k=0 



(9) 



The way to calculate coefficients afe(A) will be determined later. 

Proof. Let us divide the integral ([T]) into two integrals. Integration domain of the 
first one is the half- vicinity U of the point and integration domain of the second one is 
a remainder of integration domain of the original integral. It is easily shown by the proven 
lemma that the second integral is exponentially small as compared with exp[5'(A, a;'')]. So 
we will estimate asymptotically the first integral only. 
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In the expansion of the function S*{X,y',ip{y')) in hne with the condition ([8]) hnear 
components are absent, because the point ?/' = is a point of maximum of the function 
S*{X,y',(p{y')) in the region U'. 

Let us choose U in accordance with inequahties ip{y') < yn ^ S > at y E U*. 
Then we can represent the integral ([T]) within exponentially decreasing terms: 



F{X) = J nX,y)exp[S*{X,y)]dy. 



(10) 



Let us rewrite integral in Eq. ( ITOj) in the following way: 

F{X) = J ^X,y')dy', 

U' 



with 



exp[S*{X,y)]f*{y)dy„ 



:iii 



The integral fill I) is one-dimensional, and the function S*{X,y) reaches extremum on 
the boundary = (p{y'). Asymptotic expansion of this integral can be found through 
integration by parts. After + 1 integration we obtain the sequence: 



N 



$(A,y') = E^' 

k=0 



S*'{X,y) 
with is a unit operator and 



exp [S'*(A, 



.5 



M 



N 



S*'{X,y) 



exp[S*{X,y)]dyn 



S*'{X,y)dy^^ • 



The substitution of ?/„ = ip{y') provides the main contribution to the asymptotics, the 
value of yn = S is exponentially small as compared with the previous. Further integration 
under these conditions gives the following expansion for the function ( ITTi) : 



So 



$(A,y') = -exp[5*(A,y',^(yO)]E^' 

fc=0 



oo „ 

F{X) = - Y. J exp[S*{X,y',v{y'))]M' 



nX,y',^{y')) 



k=0 



U' 



S*'iX,y',ipiy')) 

r{X,y',v{y')) 
S*'{X,y',^{y')) 



dy' 



(12) 



The point y' = is an internal point of maximum of the integrand in the expression 
( fT2l) . Functions S*{X,y',ip{y')) and f*{X,y',ip{y')) satisfy the conditions of lemma ([21) 
and theorem ([7]), ([H]), therefore we can apply a formula for asymptotic estimation of the 
n-tuple Laplas integrals [TU] and obtain corresponding asymptotic expansion of the 
integral in Eq. (|T^ . Thus the theorem is proven. 

Note, that in the general case it is rather difficult to write a compact formula for the 
coefficients afc(A) in Eq. ([9]). They can be obtained such kind of way in any specific cases. 
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As an example, these coefficients will be obtained and written out explicitly for double 
integrals in the consideration of the asymptotic estimation of electromagnetic deuteron 
form factors in Sec. 4. Here we write out only the first asymptotic term from Eq. in 
the X variables: 

F(A) ~ -(27r)"^exp[5(A,a:°)](^^^^^^) ' \ det B\-^ /{X, x') , (13) 
where n , B are defined by conditions ^ and 

3 Electromagnetic deuteron form factors in the 

nonrelativistic and relativistic impulse approxima- 
tion 

In the nonrelativistic impulse approximation known formulas for electromagnetic deuteron 
form factors can be rewritten in the following way [H] : 



G^''m = Y. f edkk''dk'ui{k)g'^c{k,Q\k')uiik') , 
1,1' 

2 A/f ^ r 

Gq\Q') = -^Y. k'dkk''dk'ui{k)gQ{k,Q\k')u,{k'), 
^ 1,1' 

GZ''iQ') = -MdYl I k^dkk"dk'ui{k)~gl^Uk,Q\k')uiik') . (14) 
1,1' •' 

Here ui{k) are the deuteron wave functions in momentum representation, I , I' = 0,2 
are orbital angular momenta, gljiik ,k') , i = C,Q, M are nonrelativistic free two- 
particles charge, quadrupole and magnetic dipole form factors, M^^ is the deuteron mass. 
Formulas for are given in [8j. 

Let us discuss briefly possible types of the model deuteron wave functions. There are 
several classes of the deuteron wave functions: obtained with microscopic model Hamil- 
tonians of the A^A^-interaction in the non-relativistic nuclear physics (for example, see 
^J), deduced from scattering amplitudes in the Bethe-Salpeter approach and its various 
quasipotential reductions (see [12]), wave functions of the Poincare-invariant quantum 
mechanics (as an example see wave functions in the instant form of PIQM [H El El [3, [8]), 
and also wave functions calculated in the various statements of inverse scattering prob- 
lems [131 El [IS] . But independently of the method any wave function can be represented 
as the following Laguerre polynomial expansion [TSj : 



ui{k) = Yl ^'"^y r(m + / + 3/2) ' ' (^0^ ) ^ ' (1^) 



m=0 

or in the coordinate representation: 
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here Lj^^/^(x) are generalized Laguerre polynomials, T{x) is an Euler gamma function, 
the dimensional parameter tq can be related to the deuteron matter radius (see Sec. 5). 

The wave function representation as a Laguerre polynomial expansion ffT^ is very 
useful for the calculation of the asymptotic behavior of the form factors. However, one 
can avoid such representation and obtain the asymptotic expansion directly for the initial 
wave function. 

Generally, at high transferred momentum it is necessary to take into account rela- 
tivistic corrections in the electromagnetic deuteron structure. In our paper relativistic 
description of the deuteron is constructed in the framework of instant form of Poincare- 
invariant quantum mechanics (PIQM), developed by authors previously [H El El [71 E] ■ In 
this approach we present electromagnetic deuteron form factors by analogy with nonrela- 
tivistic case f|T^ . Corresponding formulas in the relativistic impulse approximation were 
obtained in our paper [B]: 

1,1' 

2 A/f ^ r 
^ 1,1' 

Gim = -M,Y.I dV~sdv^Ms)9'^Mis,Q\s')Ms'), (17) 
1,1' 

where (pi{s) are the deuteron wave functions in sense of PIQM, gQ-{{s ,Q'^ , s'), i = C,Q, M 
are relativistic free two-particles charge, quadrupole and magnetic dipole form factors. 
Formulas for free form factors are given in [8j. 

The deuteron wave functions in sense of PIQM are solutions of eigenvalue problem 
for a mass squared operator for the deuteron (see, e.g. [1]): M| {iIj) = Mj \ip). An 
eigenvalue problem for this operator is coincident with the nonrelativistic Schrodinger 
equation within a second order on deuteron binding energy e^/(4M), the value of which 
is small (M is an averaged nucleon mass). So the deuteron wave functions in sense of 
PIQM differ from nonrelativistic wave functions by conditions of normalization only. In 
the relativistic case the wave functions are normalized with relativistic density of states: 

roc rl U 

I fiik)^-7=^ = l, Mk)=V~skum, s = 4{e + M'). (18) 



1=0,2 



Nonrelativistic formulas Eq. (HM can be obtained from relativistic ones ( |T71) in the 
nonrelativistic limit. This limiting procedure can be performed in the most natural way in 
the instant form of PIQM. The reason is that in papers [H El El H E] we have constructed 
the successful formalism of the instant form of PIQM. In the case of other forms of PIQM 
(point and front forms) the obtaining of nonrelativistic limit is much more difficult. 

For obtaining the asymptotic form factors behavior at high transferred momentum in 
the nonrelativistic and relativistic cases it is necessary to estimate asymptotically double 
integrals 0141) and fll7p at Q"^ ^ oo. Notice that integrands reach its maximum value at 
the integration domain bound, and this point is not a point of extremum. In the previous 
Section the theorem defining asymptotics of n-tuple integrals of such kind was proven. 
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4 Asymptotic expansion of the deuteron form factors 



We start asymptotic expansion of the deuteron form factors from the nonrelativistic 
case. It is caused by the simphcity of the nonrelativistic formulas, so the calculation of the 
asymptotics is more clear. In what follows the relativistic calculation will be presented 
analogous to the nonrelativistic one, although more combersome. Moreover nonrelativistic 
calculation is interesting because nonrelativistic formulas for the form factors ffT^ are 
conventional, that is why its correct asymptotic expansion has universal significance. Let 
us emphasize also that the relativistic expressions for form factors and, therefore, their 
asymptotic estimations depend on the choice of the method of relativisation of the two- 
nucleon model. Nonrelativistic calculation is also of interest because it helps to clarify 
the role of relativistic effects in the electromagnetic structure of the deuteron at the 
asymptotical momentum transfers. 

As we have seen in Sec. 3, the deuteron form factors in the nonrelativistic impulse 
approximation can be represented by double integrals dHl). We will find its asymptotic 
expansion using the theorem of Section 2 and use as an example the asymptotics of the 
charge form factor. We shall estimate only the I = I' = term in the sum ffT^ because 
the asymptotics of the other terms of form factors ffT^ can be derived analogously. 

Let us rewrite the corresponding I = I' = term of the charge form factor ( fT^ using 
Eq. (USD: 

/ 9oc{k , , k') exp [S{k, k')] dk k' ^ dk' X 



We have denoted in ([19 



S{k,k') = -'^{k'' + k'''). (20) 
The expression for g^cik , , k') is commonly accepted (see, e.g., [8j): 

~g'^c{k.Q\k') = ^\e[k'- k-^)-e[k'-k-^)\{G%m + Gim). (21) 



G'e"'{Q^) are electric form factors of proton and neuteron respectively, 9{x) is a step 
function. 

In the case under consideration the space dimension n = 2, (xi, X2) = {k, k'), X = 
is a large positive parameter. Integration domain is determined by ^-functions in Eq. fl2Tl) 
and shown in Fig.l. The location of the point of maximal value of the function S can be 
obtained by analysis of ([201) and ([21]): (A;°, k 



/0\ 



454. 



Let us perform the transition to the new basis as we have descripted before. We 
perform the shift of the origin of coordinates to the point of maximal value of the function 
S. Then we rotate the obtained frame for the internal normal to the boundary in the 
new origin to be coincident with the last basis vector of the new frame. This procedure 
is illustrated in Fig.l. 

In the other words we perform the transition to the new variables in Eq. (fT9l) : 

fc = i=(t' + t) + |, A;' = i=(t'-t) + |. (22) 
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Figure 1: The integration domain, location of the point of maximal value, and transition 
to the new variables for the nonrelativistic case 



At this transformation function S{k, k') gets dependence on large parameter 



Q 



(23) 



Functions S{Q'^, t, t') , Qoci^ ; ; ^O? ^^'^ boundary of the integration domain satisfy 
the conditions of theorem 1°, 3°, 4°. The location of the point, that satisfies the conditions 
dS]), dS]), can be obtained by a simple analysis of the function f pHl) : (t°,t'°) = (0,0). Let 
us show that this point satisfies the conditions 2° of the theorem. 
It is obvious, that in the point of maximal value 



dS_ 

dn 



dS* 



dt' 



(t,i')=(o,o) 



rjQ_ 
2V2 



+ 



So the condition ^ is satisfied. 

Let us calculate now the B matrix from the condition ([U]). In our case the tan- 
gent to the domain of integration boundary vector in the point of maximal value is 
i = (l/\/2, -I/V2), i.e. the B matrix is a number: 



< 0. 



(24) 



(i,t')=(o,o) 



We note, that B is negative-definite, i.e. the point {t^,t^') = (0,0) is really the point of 
maximal value. So the point (0,0) satisfies the condition 2° of the theorem. 

So integral (fT9l) satisfies the requirements of the theorem proven in the Sec. 2. There- 
fore we can apply the asymptotic formula (IHl). 
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Calculating by analogy the other terms of the sum iHM we obtain asymptotic expan- 
sions of deuteron form factors in the nonrelativistic impulse approximation: 



„2n2 oo UNR 
"■2m 



1, 



"'2m 



k=0 'o 



k 

2m+2k+3 51 ^fcP 
p=0 



'9 /9\P+''+^ fi2m 



Q 



dt 



2m •' * 



(25) 
(26) 



i=0 



ho = 1, hp = h~ip- {k + p- l)h-ip-i, hk = (-i)'''(2A; - 1)!! 
frit,Q',t') = A, E k^'-'k^^'udk) &,Q',t')u,{k') , 



Ll'=0,2 



with k = k{t,Q^,t') ,k' = k'{t,Q^,t'), variables t, t' are denoted in ([221), i = C,Q,M, 
Ac = 1, Aq = 2MJ/Q2 , Am = -M,, 



ui II is defined by equalities: 



r"^'"^(t,Q^^') = ^/r(t,Q^t'). 



2, 2 



'Uo(^) = Mo(^) e 2 ^ ■U2(^) = M2(^) A; e 2 . 



(27) 



Let us perform now the calculation of the relativistic asymptotics of deuteron form factors. 
To estimate asymptotically integrals (fT7j) we proceed analogously to nonrelativistic case, 
i.e. we use relativistic analogs of corresponding nonrelativistic formulas ( fT9l) - (!Ml) . Now 
the free relativistic charge form factor in ( JT71) at / = /' = is given in Ref. [8]: 

gllis, Q\ s') = Ris, Q\ s') Q'[{s + s' + Q') {G^Q') + GUQ')) 



f28) 



^^E^AiiQ"^) electric and magnetic form factors of proton and neutron respectively. 



9(JE = ( o "^os iUi cos UJ2 + - sin Ui sin UJ2 I , g'cM = ( o ^'^^ sin ^2 — - sin cji cos Ci;2 



6 



00 



1 



6 



R{s,Q^s'' 



{s + s' + 



Ui and UJ2 are angles of the Wigner spin rotation. 



s - 4M2)(s' - 4M2) [A(s, 5')]^^' y^l + QV4M2 
e(s, g^ s') = Jss'Q^~M^\{s,-Q\s') , 



Wi = arctan ■ 



UJ2 = arctan 



a(s,s')e(s,Q2,s') 



M(s + s' + Q'^)a{s, s') + Vss'(4M2 + q2 



(29) 
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where a{s, s') = 2M + y/s + V^, ^{s, Q^, s') = 6{s' — si) — 6{s' — S2), 6' is a step function, 
A(a, 6, c) = + 6^ + - 2{ab + ac + be), 



«i,2 = 2M^ + ^(2M^ + g2)(. - 2M2) T ^\/^'(^' + " 4^') • 

To obtain relativistic asymptotic expansion we also perform transition to the new basis 
(shift and rotation). The function S and the boundary of the integration domain differ 
from nonrelativistic ones, so it is necessary to perform a special analysis. In other words, 
instead of change of variables (1221) we perform the following replacement: 



s = -^\^t' + 2M^ + M^Q^ + 4M2 , s'=-^ \^t' - g j + + M^Q^ + . 

(30) 

Then we obtain asymptotic expansion of the relativistic deuteron form factors by 
analogy with nonrelativistic case: 

,,2 CO 1 

Gf (Q2) ^ ^-^iM^Q^+4M^-2M^) ^ p^T{m + -), (31) 

=0 (2^^)- 2 



m= 



oo 2m / -I \p r,|A:+4m-2p+| Q2m— 2p 



t=0 

(32) 



/f(t,Q^t') = ^. E u,{k)glM\t')u,{k') ^''^ M2)^^(^/4 M2) ^ ^ 

/f('")(t,Q^t') = ^/f(t,Q^t')• 

Functions = /c(s) , k' = k'{s') are specified in Eq. ffTHj) . s = s(t, Q^, t') , s' = s'{t, Q'^, t'), 
variables t, t' are denoted in flHUj) . C|™ are the binomial coefficients. 

Asymptotic expansions ( !25|) and ( I3T1) are convergent power series in inverse degrees 
of the parameter Q with known coefficients. The asymptotic expansion of this type is 
obtained in this work for the first time. 

One can see from formulas (1251) and (I3ip . that relativistic corrections change the 
behavior of form factors at high momentum transfer. In particular, exponential multiplier 
index is in the nonrelativistic case, but in the relativistic case it is Q at — >■ oo. It 
seems to be a general feature of our relativistic approach to the description of composite 
systems, in particular, we have obtained the similar result in consideration of asymptotic 
behavior of the pion form factor in the composite quark model |16j . 

5 Asymptotics of the form factors for the conven- 
tional wave functions representation 

In this Section we represent the obtained asymptotical expansions ( l25l) and ( |3T1) in 
terms of initial wave functions in the left side of Eqs. (JT5l) . (1181) . For this representation it 
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is necessary to replace functions ui{k) by functions ui{k) in ( l25l) and (I3T1) using (|T8l) . (1271) . 
Keeping the main term on 1/Q in asymptotic expansions and flHTj) one can obtain 
the next asymptotic formulas in terms of functions Ui{k) and ipi{s) from ffT^ . ffTHl) : 



(33) 



t=0 
t'=0 



Ll'=0,2 



4/ I 

v s s' 



, i = C,Q,M . (34) 



t=0 
t'=0 



Let us note, that similar asymptotic representation can be obtained for any finite number 
of terms in asymptotic expansions ([2SD, (IHT]) . 

In the modern calculations the deuteron wave functions are usually represented as a 
discrete superposition of Yukawa- type terms (see, e.g., [9]): 



? Mk) = \-Yl 



^ (P + m^) 
or in the coordinate representation: 

«o('^) = E C'j-exp {-rrij r) 



TT ^ (F + mj) ' 



(35) 



'^2{'r) = E -DjCxp {—rrij r) 



3 3 
1 + + 



rrij r {rrij r)^ 



(36) 



nij = a + mo (j — 1) , a = \J M 



Coefficients Cj, Dj, maximal value of the index j and mo are determined by the best fit 
of corresponding solution of Schrodinger equation. 

The deuteron wave function analytical form ( l36l) results in the right behavior of the 
wave functions at large distances: 



Uo{r) ~ exp(— ar), U2{r) ~ exp(— ar) 1 



3 3 

+ 



(37) 



(38) 



(39) 



Let us substitute the wave functions (l35l) to ( l33l) and ( IMI) . and then obtain the first 
asymptotic terms of the nonrelativistic deuteron form factors: 



The deuteron wave functions behavior at small distances: 

■Uo(r) ~ r , U2{r) ~ , 
is provided by imposing the following conditions on coefficients Cj and Dj: 



D 



^ = 

2 • 
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1 2 



16 



NR 



EC. 



NR 



3M, 



1 2f 



4 



1 2^^Md 



EC: 



G',,m + Cl,{Q' 



(40) 
(41) 
(42) 



The dimensional parameter ro can be found from the expression for the deuteron 
matter radius in our deuteron model: 



4 Jo 



(43) 



One can substitute wave functions of the form (fT6|) into this expression. So formula (l43l) 
specifies an algebraic equation for vq. Solution of this equation can be found numerically. 

It should be pointed out that main terms of expansion of charge and magnetic form 
factors in ( HOl) - (H2l) are determined by S'-state of deuteron only. The D-wave function 
gives the contribution to the main term of the quadrupole from factor. Its faster decrease 
at ^ oo in comparison to the other form factors is a consequence of a faster decrease 
of a D-wave function at small distances in comparison to 5*- wave f l38l) . From the mathe- 
matical point of view the type of leading terms in fHni) - fH^ is a consequence of conditions 
on the coefficients ( l39l) . The modification of these conditions obviously results in change 
of the main terms in ( l40i) -( H2l) . From these formulas it is also noticed that asymptotic ex- 
pansions for the deuteron form factors contain dependence on the asymptotics of nucleon 
form factors. 

We emphasize, that in the other deuteron asymptotics investigations only the power 
dependence on the transferred momentum was calculated as a rule. In the present paper 
we give a rigorous calculation of a multiplicative preasymptotical constant. 

One can calculate relativistic asymptotics of form factors by analogy with nonrela- 
tivistic case. For this calculation we use the formulas (fT^ . f lMl) . flH^ . As a result we 
obtain: 



Cc,m{Q^ 



'Q 



22M3 



(44) 



(45) 



Notice that asymptotic expansions 



2-M4 

and fH^ . (H31) are obtained for the first time 
in our work. It is interesting to compare obtained asymptotic estimations fHUj) - fH^ . 
f H4l) .f H5|) with observable behavior of the deuteron characteristics. At present time there 
exists the experimental information about function AIQ"^) entered the differential cross 
section of the elastic erf-scattering. This function is expressed in terms of the deuteron 
form factors pQ. The values of function A(Q^) are known up to — 6 (GeV/c)^. For 
the comparison with experimental data one needs to specify asymptotics of the nucleon 
form factors. It is naturally to choose for nucleon form factors the asymptotic which 
is predicted by the quark model G^" ~ 1/Q^- Under these conditions the power 
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dependence on of the function A{Q'^) coincides with experimentally observed one. 
The physical consequences will be examined in detail in the other paper. 

6 Conclusion 

The theorem defining asymptotics of multiple integrals of some special type is proved. 
With help of the proven theorem the asymptotic expansion of the deuteron electromag- 
netic form factors at oo is calculated for the first time. The expansion is repre- 
sented as a convergent series on in inverse powers of momentum transfer. The asymptotic 
of the form factors is found in terms of the conventional representation of the deuteron 
wave function as a discrete superposition of Yukawa-type terms. The asymptotic behavior 
of the form factors is calculated in the nonrelativistic impulse approximation and in the 
relativistic invariant impulse approximation proposed by the authors in the instant form 
of the Poincare-invariant quantum mechanics previously. It is established that relativistic 
corrections change the power dependence of the form factors on the momentum transfer 
sX oo and slow down its decrease. It is also found that relativistic effects result in 
the agreement of the theoretical asymptotics and the experimentally observed behavior 
of the structure function A{Q'^) at highest achieved momentum transfers. 
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